Attractions between charged colloidal spheres mediated by correlated absorbed ions 



O 

o 
o 

> 

O 



■4— » 
•4-* 



o 
o 



> 
in 

O 
O 

o 
o 

-I— > 



T3 

c 

o 
o 



x 



Francisco J. Solis 1 and Monica Olvera de la Cruz 2 
1 Physics Department, Hiram College, Hiram OH, 44238 
2 Department of Materials Science and Engineering, Northwestern University, Evanston IL 60208 
e-mail: 1 solisfj@hiram.edu, 2 m-olvera@northwestern.edu 

We describe arrangements of ions capable of producing short-range attractive interactions be- 
tween pairs of charged colloidal spheres in the low temperature strongly correlated limit. For parti- 
cles of radius R with bare charge Z and comparable absorbed charge — N (N ~ Z), the correlations 
contribution to the spheres self-energy scales as N 3 ^ 2 /R, and as N/R for the interaction energy 
between two touching spheres. We show that the re-arrangement of charges due to polarization 
plays an insignificant role in the nature and magnitude of the interaction. 



Electrostatic attractions between colloids have been re- 
ported in a large number of circumstances [jj]. Experi- 
ments in colloidal suspension confined into a small (al- 
most two dimensional) region |2|-||] have shown direct ev- 
idence of these attractions. Simulations of colloidal sus- 
pensions in bulk and confined systems J?]] have also 
reported electrostatic attractions. The simulations have 
explored, in particular, the strongly correlated limit in 
which the thermal energy of the system is smaller than 
the typical interaction energy between the charges ab- 
sorbed to the colloids. 

Recent studies on colloidal suspensions have elucidated 
the effect of depletion forces and of Van der Waals 
interactions pL as well as the importance of including 
correlations to predict attractions. In this let- 

ter we derive analytic expressions for the interaction en- 
ergy between two colloids in the strongly correlated limit 
by considering the distribution of absorbed counterions 
around them as they approach. Though correlations be- 
tween absorbed counterions are known to induced short 
range attractions between charged surfaces [ |F2|jl3]] , and 
between aligned rods Jl^-[l^1, it has been noted || that 
the geometry change to colloids is not trivial. 

Systems involving large charged objects (polyions) 
whose charge is compensated by mobile low-valence ions 
can be study by considering first the electrostatic ground 
state of a subsystem that consists of one or few polyions 
with a fixed number of counterions. Once the ground 
state is constructed, is possible to draw conclusions about 
the system by introducing thermal fluctuations and set- 
ting the subsystem in contact with the rest of the system 
which then acts as a thermal and particle reservoir. This 
approach is justified by the fact that there is a small re- 
gion surrounding the polyion in which the electrostatic 
interactions are so dominant that thermal fluctuations 
are only a minor correction. It has been used, for exam- 
ple, to study the interactions between charged surfaces 
p3[ and between charged linear polyelectrolytes p3 , n5| , 
the conditions of solubility of DNA and other linear poly- 
electrolytes in multivalent salt solutions |l7|]Ts|], and to 



describe the differences between charging colloids or sur- 
faces with monovalent or multivalent counterions [ jl9| . 
Results from these theories were used in Ref. || to ob- 
tain simple but accurate predictions that explain many 
colloidal simulations results. 

While in principle it might be possible to solve the 
ground state of a system consisting of one charged sphere 
and a set of N absorbed point-charges, such a ground 
state has to be solved almost independently for each dif- 
ferent N. It would be even more difficult to obtain exact 
results for the case of two spheres, although numerical 
results can be obtained for N not too large. Instead of 
attempting such solutions, we propose an approximate 
solution that satisfies a number of desirable physical con- 
ditions. We first decompose the charge distribution on 
the sphere into a multipolar series for each sphere. We 
then truncate the series and propose values for the dom- 
inant terms. Using this ansatz, we obtain an approxima- 
tion for the interaction energy between two spheres. 

Our main results are simple to understand. We recover 
the expected negative self-energy for a single charged 
sphere due to the discrete nature of the absorbed coun- 
terions, which is of the form ~N 3 / 2 /R M. Indeed, since 
R ~ iV 1 / 2 a, where a is the typical distance between the 
absorbed charges, this self-energy is roughly proportional 
to —N/a, as in the case of a flat surface. If the spheres 
have bare charge Z, and N absorbed (opposite) charges, 
the effective charge is Q = Z — N, and thus there is an 
interaction energy between two spheres at a distance r 
of the form Q 2 /r, but there are important contributions 
arising from the particle nature of the absorbed charge. 
When the spheres are in contact there is a region in which 
the charges from both spheres are strongly correlated. 
In this region the bare and absorbed charge density are 
doubled and the typical distance between charges is now 
of the order of <z/v2- The charges in both spheres re- 
main correlated up to when the separation between cor- 
responding points in the spheres are separated by a dis- 
tance of the order a. This defines an effective area of 
contact between the spheres given by ir Ra, that contains 
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Na/4R charges. The interaction energy is then negative 
and of order N/R. Thus, there can be a very large attrac- 
tive component to the interaction between the spheres 
capable of producing a net attraction between them for 
almost neutral spheres (N ~ Z). 

We assume that the total bare charge of the colloidal 
sphere Z is uniformly distributed on its surface or, equiv- 
alently, concentrated at the center of the sphere. The 
surface charge density of the N absorbed particles is then 
written as p(r) — £j£(r — rj) where Tj are vectors point- 
ing from the center of sphere to a position at the surface. 
The multipolar moments M™ for this distribution are 



M m = 
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where Y™ are spherical harmonics. The charge density 
can be then re-expressed as the multipolar expansion: 
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Using Maxwell's representation |g0| of the spherical har- 
monics it is also possible to replace the surface density by 
a singular distribution placed at the center of the sphere 
given by 



P(r) = £ 
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We will use the more compact notation p = 
^M™{n, to} where the brace {n,m} stands in for the 
suitable functional form in this expression. The potential 
generated by this charge density outside of the sphere, at 
r > R, is 
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In the low temperature limit we expect the charges 
to occupy positions that up to permutations and rota- 
tions of the whole sphere are uniquely determined. For 
large N, the inhomogeneities arising from the particle 
nature of the charges in this static ground-state distri- 
bution should create only negligible contributions for the 
first few multipolar moments other than n — 0, with 
the next non-vanishing moments arising at n values that 
start to resolve distances of order a, namely n ~ TV 1 / 2 . 
To make analytical calculations feasible we truncate the 
multipolar expansion and represent the sum of the N 
delta functions as a smooth function with N clearly dis- 
tinguished maxima. Since the interaction between the 
absorbed charges is repulsive one can expect four nearest 
neighbors per particle in the ground state and we also 
require this from our truncated expansion. Our ansatz is 
then a truncation of the series Eq.(^) that only consider 
moments with n values up to n = 31/2 with I = \/2N. 



Further, within this range we take only four moments to 
be non zero: 

p = -N{0, 0} + V I.,, {I, 0} + 
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The precise values of the moments given here are dis- 
cussed below. 

The proposed charge density Eq.(||) has N maxima, 
distributed between I = \J N/2 meridians generated by 
the to = I terms. Each meridian has 1/2, non uniformly 
distributed maxima generated by the to = term. A 
consistent extension of our ansatz would give the next 
non-zero moments at n and ±to values that are integer 
multiples of Z. In the limit of infinite radius and near the 
equator of the sphere the distribution becomes equiva- 
lent with the description, by a few Fourier modes, of a 
square lattice of charges on a plane. This limit is indeed 
the inspiration for our choice. The major inaccuracies 
of this distribution are near the poles where the maxima 
are closer than they should be. 

For large N, and near the equator, the spherical har- 
monic YJ° has an oscillatory form tt^ 1 cos(W + ip), with 
9 the polar coordinate and ip a phase. If we make the 
position of the charges to coincide with the peaks of 
this approximate distribution, we obtain, using Eq.(Q), 
M ; ° = n- 1 / 2 (2N) 3 / 4 R l . Again, near the equator, we 
have that Y/^ 1 = e ±d ^f(9), where <fi is the azimuthal an- 
gle, and / is a multiple of the Legendre function P\ that 
has a maximum in the equator and is important only in 
a narrow region. We assume that the particles positions 
coincide with the maxima of the real parts of the waves in 
the (f> coordinate while the sum over their contributions 
along the is approximated by a multiple of the integral 
/ f{9)d9 from which M\ = 2 7 / s (N) 5 / s R l /(tt 3 / 4 ). 

To discuss the interaction between spheres we must 
consider the redistribution of charges that occurs when 
they are put in close contact with each other. To first 
order we can represent this redistribution as a flow gen- 
erated by the projection, onto the sphere, of a uniform 
displacement of charges in the direction of the polar- 
ization. If the displacement of the charges is d, the 
new charge distribution is approximately p(r) + Ap(r) = 
p(r) + d- Vp(r). Using the representation Eq. ^ and as- 
suming the displacement to be in the x direction, it can 
be shown that the displacement applied to our ansatz 
Eq.(||) gives rise to new non-zero extra moments: 
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This equation displays only the positive m terms. The 
magnitude of the generated dipole moment is P = dN, 
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but as this expression shows, the flow of charges also 
affects the higher order multipoles. The linear approxi- 
mation is valid only if the average displacements of the 
charges are smaller than the typical distance between 
them d < a 

The energy of interaction between the absorbed 
charges and the colloid is —ZN/R. If the decomposi- 
tion of the ground state configuration of the absorbed 
charges into multipolar moments is known, their self 
energy is E = J2 n ,m M n M n m /2R 2n+1 . This expression 
is infinite since it contains the self-energy of every ab- 
sorbed charge. The harmonic decomposition of a single 
point charge shows that its self-energy splits into equal 
contributions of magnitude 1/2R for each n number giv- 
ing, symbolically, E sing i e = (l/2R)J2 n 1 - Subtracting 
TV times this amount from the energy at each n level we 
obtain: 
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Evaluating the energy of our ansatz for the charge dis- 
tribution and including the possibility of a charge redis- 
tribution with dipole moment P the main terms in the 
evaluation of the energy are 



Z 2 -2ZN + N(N - 1) 
2R 
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where a = (3/2 3 / 2 ) - 1/(tt2 1 / 2 ), and (i = 1/{2 1 / 2 tt). In 
the second and third terms we have retained only the 
contributions with highest powers in N . The monopole 
contributions have been explicitly separated. In the ab- 
sence of polarization, this result is clearly consistent with 
the expected behavior of the self-energy in the N — > oo 
limit that is equivalent to charges in a flat surface in 
which the self-energy scales as N(l/a) ~ N(N/ R 2 ) 1 / 2 . 

To consider two spheres we orient the reference systems 
for each of them so that their centers lie in the x-axis of 
both systems. The total distance between centers is r. 
The interaction energy can be written as the integral of 
the potential generated by one sphere <f> times the charge 
density of the second sphere p' , E — J dr<j>p' . Writing 
down both <f> and p' in their multipolar expansions and 
using the Maxwell representation, the integral becomes 
a sum of terms that are products of derivatives of 1/r 
times derivatives of a delta function at the center of the 
second sphere. Repeated integration by parts leads to 
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where the coefficient is, using s = n + m,t = n — m, s' — 
p + q,t' = p — q, S = s + s' ,T = t + t' . 
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The point of evaluation of the spherical harmonic reflects 
the direction of the relative position of the spheres in 
the chosen coordinates, namely the x axis. It can be 
shown that for large moment numbers the magnitude 
of the coefficient has its maximum values when n = p 
and m — ±q, and decays away from this situation as 
exp(-(s - (S/2)) 2 - (t - (T/2)) 2 ). Therefore, the main 
contributions to the energy from high order multipoles 
appears when the multipoles in both spheres correspond, 
that is, when the charge distributions are similar and in- 
terlock. In our truncation we are able to pick the relative 
phases of the multipoles in the spheres to get an attrac- 
tive energy. 

We will consider only the case of two spheres of equal 
radius. Since we are interested in exploring the possibil- 
ity of attractive interactions between spheres we consider 
only the case in which they touch r = 2R. If attractions 
are present this will be their optimal position. Consider- 
ing the case of spheres that have equal total charge and 
magnitude of polarization we obtain the following inter- 
action energy: 

_{Z-N) 2 N (P-P')(Z-N) 
mt ~ 2R 7 2i? {2Rf + 
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energy is obtained 



from adding the self-energy and the interaction energy. 
Minimizing the total energy with respect to the polar- 
ization we arrive to the conclusion that the polarization 
of the spheres should be antiparallel P = — P', and with 
magnitude P — RQ/ (3N 1 / 2 . The polarization is negligi- 
ble in the limit of large N. The reason for this is that 
a flow of mobile charges towards the point of contact re- 
quires a large increase of energy in a large depleted area, 
while the decrease in energy occurs only in the small re- 
gion of contact. Thus, we find that the ground state is 
very rigid, and the interaction between spheres does not 
alter much its local structure. The polarization correc- 
tions can then be ignored and the interaction energy is 



simply the result given in Eq.(ll) without polarization 
terms. If the effective charge is small, the interaction 
is dominated by the high multipole correlations and it 
scales as N/ R. The rigidity of the local structure makes 
this result robust when applied to considerations of more 
than two spheres. 

Given our result for the evaluation of the interaction 
energy, we are ready to comment on the problem of at- 
traction in more realistic conditions. The Bjerrum num- 
ber B = e 2 /aAireoksT compares the electrostatic inter- 
action between charges at a distance a with the typical 
thermal energy ksT. We can consider z-valent counteri- 
ons and then, if we have z 2 B > 1, the correlation energy 
can be very strong. The exact number of charges in the 
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spheres is controlled by the equilibrium of the energy be- 
tween absorbed charges and a chemical potential /i for 
the extraction of them from the volume not occupied by 
the colloidal spheres. To solve for the number of absorbed 
charges we can minimize the free energy F — E(N) — jiN 
with respect to N. Short range attractions between two 
charged colloids proportional to N/R result from the "in- 
terlocking" correlations of the thin layer of counterions at 
the surface of the colloids in the area of contact. These 
correlations easily arise, for high Bz 2 , when the num- 
ber of absorbed counterions is large and comparable to 
the bare charge. The relative location of spheres influ- 
ences the amounts of absorbed charges but their specific 
values are mostly dependent on the interaction parame- 
ter Bz 2 and on the chemical potential /i. This suggests 
that it is more likely to observe attractions in solutions of 
high valence counterions or in systems where the acces- 
sible volume for the counterions is small as in the case of 
colloids in confined geometries. Though in principle the 
confinement of colloidal suspensions by weakly charged 
surfaces imposes non-trivial boundary conditions that af- 
fect the local electrostatic fields, we can always analyze 
the local configurations of ions under these externally im- 
posed conditions. If the external conditions are successful 
in creating attractions, the local arrangement of charges 
are as we have described here. 

The explanation for the presence of long-range attrac- 
tions put forward in Ref. j|] is that there are attractive 
metastable states characterized by pairs of spheres with 
opposite total charges. The overcharging of some spheres 
is possible due to the negative correlation self-energy in 
each sphere Eq.(^). What our calculation shows is that 
even if the effective charges are equal, the ground state is 
attractive. The metastable states, on the other hand, re- 
duce the potential barrier to achieve close contacts, and 
at the free energy level might indeed provide a long range 
attractive interaction. At short distances we expect the 
charge to redistribute to the type of configurations dis- 
cussed here. 
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